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Abstract 

The EGRET telescope has identified a gamma-ray source at the Galactic center. 
We point out here that the spectral features of this source are compatible with the 
gamma-ray flux induced by pair annihilations of dark matter weakly interacting 
massive particles (WIMPs). We show that the discrimination between this interpre- 
tation and other viable explanations will be possible with GLAST, the next major 
gamma-ray telescope in space, on the basis of both the spectral and the angular 
signature of the WIMP-induced component. If, on the other hand, the data will 
point to an alternative explanation, we prove that there will still be the possibility 
for GLAST to single out a weaker dark matter source at the Galactic center. The 
potential of GLAST has been explored both in the context of a generic simplified 
toy-model for WIMP dark matter, and in a more specific setup, the case of dark 
matter neutralinos in the minimal supergravity framework. In the latter, we find 
that even in the case of moderate dark matter densities in the Galactic center region, 
there are portions of the parameter space which will be probed by GLAST. 
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1 Introduction 



Experimental cosmology has been steadily progressing over the latest years. 
The emerging picture has been recently reinforced by the data from the Wilkin- 
son Microwave Anisotropy Probe (WMAP) which have pinned down several 
fundamental parameters to a remarkable level of precision. In particular, in the 
latest global fit [1] , the contribution to the critical density of non-relativistic 
matter has been found in the range VL m h 2 = 0.135t(j;oo8 (here h is the Hubble 
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constant in units of 100 km s 1 Mpc h — O.Tllo^fl]), m uch larger than 
the baryonic term, Vt b h 2 = 0.0224 ± 0.0009. 

Unveiling the nature of non-baryonic cold dark matter (CDM) is one of the ma- 
jor challenges in science today. Weakly interacting massive particles (WIMPs) 
are among the leading dark matter candidates. They would naturally appear 
as another of the thermal leftovers from the early Universe, and, at the same 
time, their existence is predicted in several classes of extensions of the Stan- 
dard Model of particle physics. The most popular of such candidates is the 
lightest neutralino in R-parity conserving supersymmetric models. Consider- 
able effort has been put in the search for dark matter WIMPs in the last 
decade, with several complementary techniques applied (for a recent review, 
see, e.g., [2]). Among them, indirect detection through the identification of the 
yields of WIMP pair annihilations in dark matter structures [3, 4] seems to be 
a very promising method. In particular, we will focus here, as a signature to 
identify dark matter, on the possible distortion of the spectrum of the diffuse 
7-ray flux in the Galaxy due to a WIMP-induced component, extending up 
to an energy equal to the WIMP mass (a list of other recent analysis on this 
topic includes Refs. [5,6,7,8,9,10,11]). 

The EGRET telescope on board of the Compton Gamma-Ray Observatory 
has mapped the 7-ray sky up to an energy of about 20 GeV over a period 
of 5 years. The data collected by EGRET toward the Galactic center (GC) 
region show [12] high statistical evidence for a gamma-ray source, possibly 
diffuse rather than point-like, located within 1.5° of the GC (/ = b = 0°). The 
detected flux largely exceeds the diffuse 7-ray component expected in the GC 
direction with a standard modeling of the interaction of primary cosmic rays 
with the interstellar medium (see, e.g., [13]); the latter fails also to reproduce 
the spectral shape of the GC source. Although other plausible explanations 
have been formulated, it is very intriguing that the EGRET GeV excess shows, 
as basic features, the kind of distortion of the diffuse 7-ray spectrum one would 
expect from a WIMP-induced component, assuming that the dark matter 
halo profile is peaked toward the GC. We will identify for which classes of 
WIMP compositions and masses, fair agreement with the measured flux can 
be obtained. 

No firm conclusion about the nature of the GC excess can be driven from 
data available at present; on the other hand, the picture is going to become 
much clearer in the near future. The Gamma-ray Large Area Space Telescope 
(GLAST)[14] has been selected by NASA as the next major 7-ray mission, 1 



A list of people and institutions involved in the collaboration together with the 
on-line status of the project is available at http://www-glast.stanford.edu For a 
detailed description of the apparatus see [15]; a discussion of the main scientific 
goals can be found in [16]. 
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and is scheduled for launch in the first half of 2006. Compared to EGRET, 
GLAST will have a much larger effective area, better energy and angular 
resolutions, as well as it will cover a much wider energy range. GLAST will 
perform an all-sky survey of 7-ray sources, with scientific objectives including 
the study of blazars, 7-ray bursts, supernova remnants, pulsars, the diffuse ra- 
diation in the Galaxy, and unidentified high-energy sources. The identification 
of dark matter sources has been indicated as one of its main scientific goals. 
We illustrate here the conditions under which it may be feasible that GLAST 
will single out a dark matter source located at the GC. 

The paper is organized as follows. In Section 2 we show that the GC EGRET 
7-ray excess can be modeled with a component due to WIMP annihilations. 
In Section 3 we discuss what kind of information GLAST could provide to 
confirm such hypothesis. In Section 4 we illustrate the perspectives of dark 
matter detection with GLAST in the case in which the EGRET excess will 
be found to be due to another kind of source. In the first part of the paper 
we discuss features for the dark matter signal in the contest of a generic toy- 
model; in Section 5 we will focus on neutralino dark matter candidates in 
the minimal supergravity framework, applying to this specific case the tools 
developed in the other Sections. Conclusions are in Section 6. 



2 A Dark Matter Source at the Galactic Center? 



In Table 1 we report the flux per energy bin for the GC gamma-ray source 
as measured by EGRET, together with the expected flux from cosmic ray 
interactions in a standard scenario [12,17] (see also Table 2 and Fig. 4 in [12]). 
As already mentioned there is a significant mismatch between the two: we 
entertain here the possibility that the bulk of the high energy flux is due to 
pair annihilations of dark matter WIMP's in the GC region. Hence, we assume 
that the total flux measured by EGRET can be described as the superposition 
of two contributions: 

• i) a component due the interaction of primary cosmic rays with the interstel- 
lar medium, with spectral shape defined by the function Sb(Ej) (background 
contribution) 

• ii) a component due to WIMP annihilation in the dark matter halo, whose 
energy spectrum is defined by S X (E^) (signal contribution). 

We write the total 7-ray flux as: 

7 = <t>b + <t> x = N bS b + N X S X , (1) 
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Table 1 

Estimated values for the Galactic diffuse 7— ray component (second column) and 
EGRET data from a region of 1.5° around the GC (third column), extracted from 
[12]. 



Energy Bin 
(GeV) 


Expected Diffuse 7— Ray Flux 
(cm-ViGeV^sr- 1 ) 


Total 7-Ray Flux 
(cm^s^GeV^sr- 1 


0.03 - 0.05 


3.7 


• 10~ 3 


(5.0 ±0.8) • 10~ 2 


0.05 - 0.07 


1.8 


• 10~ 3 


(1.3 ±0.2) • 10~ 2 


0.07-0.1 


1.1 


• 10" 3 


(6.1 ±0.5) • 10" 3 


0.1 - 0.15 


6.2 


• 10- 4 


(4.4 ±0.2) • 10~ 3 


0.15 - 0.3 


2.6 


• 10~ 4 


(2.03 ± 0.06) • 10" 3 


0.3-0.5 


1.0 


• 10- 4 


(9.5 ±0.2) • 10~ 4 


0.5-1 


3.5 


• 10~ 5 


(3.9 ±0.1) • 10~ 4 


1 - 2 


9.1 


• 10~ 6 


(1.52 ±0.03) • 10~ 4 


2-4 


2.0 


• 10~ 6 


(3.2 ±0.1) • 10- 5 


4-10 


2.3 


•10- 7 


(3.1 ±0.2) • 10~ 6 



where N b and N x are dimensionless normalization parameters, respectively, 
for the standard and exotic flux, which we define below. 

2.1 The background component 

There are three mechanisms which give rise to diffuse 7-ray radiation in 
the Galaxy: production and decay of 7r°s, inverse Compton scattering and 
bremsstrahlung (see, e.g., [13]). According to standard scenarios, in the en- 
ergy range we will consider, E 1 > 1 Gev, the dominant background source 
is given by 7r° decays. Photons are generated in the interaction of primary 
cosmic rays with the interstellar medium via: 

p + X -> . . . -> 7T° -> 2 7 

He + X -> . . . -> 7T° -> 2 7 , 

where X stands for an interstellar atom, mainly H and He. We have simulated 
the induced 7— ray yield according to standard treatments (see, e.g., [18, 19]) 
and as implemented in the Galprop software package [13]. We assume that the 
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p and He cosmic ray fluxes in the Galaxy have the same energy spectra and 
relative normalization as those measured in the local neighborhood [20,21], and 
that the He fraction of the interstellar gas is 7%. We write the background 
flux, splitting it into two factors: 

S b (E,) = 1 • Em(E 1 ) (2) 
(1 cnrsr) 



and 



b (lcm-V- 1 ) \J s 4tt 0r m (/ = O) ' 1 ' 



Here Em(E J ) [GeV -1 s _1 ] is the local emissivity per hydrogen atom, i.e. the 
number of secondary photons with energy in the range (E y , E^ + dE^) emitted 
per unit time per target hydrogen atom, for an incident flux of protons and 
helium nuclei equal to the locally measured primary proton and helium fluxes. 
The factor N b is instead associated to the interstellar hydrogen column density 
nn{l), integrated along the line of sight and weighted over the proton primary 
flux at the location /, 0£" m (Z), normalized to the local value <^" m (7 = 0). 

Above an energy of about 1 GeV the background spectrum (and therefore the 
function 4>b) can be described by a power law of the form (pback oc E~ a with 
the same spectral index as the dominant primary component, i.e. the proton 
spectral index a = 2.72. 

The relative normalization of the primary components in different places in 
the Galaxy can be estimated once a radial distribution of primary sources 
is chosen (following, for instance, the radial distribution of supernovae) and 
then by propagating the injected fluxes with an appropriate transport equation 
(this is what is done, e.g., in the Galprop code [13]). On the other hand, the 
hydrogen column density toward the Galactic center is very uncertain; we 
choose therefore to define the spectral shape of the background through the 
function Sb, and to keep TVj, as a free normalization parameter. 



2.2 The signal component 



The production of 7-rays in a dark matter halos made of WIMP's follows 
essentially by the definition of WIMP, regardless of any specific scenario one 
has in mind. The signal scales linearly with the pair annihilation rate in the 
limit of non-relativistic particles. 

We consider a generic framework in which the dark matter in the Galactic halo 
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is made of particles x> WIMP dark matter candidates with mass M x and total 
pair annihilation rate into lighter Standard Model particles av (in the limit of 
zero relative velocity). Among the kinematically-allowed tree-level final states, 
the leading channels are often bb, ti, t + t~, W + W~, Z°Z°. This is the case, e.g., 
for neutralinos and, more generically, for any Majorana fermion WIMP, as for 
such particles the S-wave annihilation rate into the light fermion species is sup- 
pressed by the factor m'j/M 2 , where m/ is the mass of the fermion in the final 
state. The fragmentation and/or the decay of the tree-level annihilation states 
gives rise to photons. Again the dominant intermediate step is the generation 
of neutral pions and their decay into 27. The simulation of the photon yield is 
standard; we take advantage of a simulation performed with the Lund Monte 
Carlo program Pythia 6.202 [23] implemented in the DarkSUSY package [24]. 

Suppose that the dark matter halo is roughly spherical and consider the in- 
duced 7-ray flux in the direction that forms an angle ip with the direction of 
the Galactic center; the WIMP induced photon flux is the sum of the contri- 
butions along the line of sight (l.o.s): 

J l.o.s x 



where Bf is the branching ratio into the tree-level annihilation final state 
/, while dNf/dE is the relative differential photon yield. The WIMP mass 
density along the line of sight, p(l), enters critically in the prediction for the 
flux, as the number of WIMP pairs is equal to 1/2 p(l) 2 j 'M£. It is then useful 
to factorize the flux in Eq.(4) into two pieces, one depending only on the 
the particle physics setup, i.e. on the cross section, the branching ratios and 
the WIMP mass, and the other depending on the WIMP distribution in the 
galactic halo. We rewrite Eq.(4) as [25]: 



W) = 3.74.10 -°( 10 - 26cnA -, ) (— j E^/ 

•J($) cm-V^eV- 1 ^- 1 (5) 

where we introduced the dimensionless function J, containing the dependence 
on the halo density profile, 

^ ( o.3 GeL-O M^'W < 6 > 



More precisely, given a detector with angular acceptance Af2, we have to 
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Fig. 1. In the left panel: differential 7-ray yield per annihilation (see Eq.(4)) for a 
fixed annihilation channel (bb) and for a few sample values of WIMP masses. For 
comparison we also show the emissivity, with an arbitrarily rescaled normalization, 
from the interaction of primaries with the interstellar medium. In the right panel: 
differential yield per annihilation for a few sample annihilation channels and a fixed 
WIMP mass (200 GeV). The solid lines are the total yields, while the dashed lines 
are components not due to ir° decays. 

consider the average of J(ip) over the solid angle AQ around the direction ip: 



An 



AQ 



(7) 



To compare with the GC EGRET gamma-ray source, we will consider AVt ~ 
10~ 3 sr, i.e. the same magnitude as the angular region probed by the EGRET 
telescope. 

As for the background component, we split the signal into a term which fixes 
the spectral shape of the flux, plus a normalization factor. In the notation 
introduced in Eq.(l), we label N x = (J('0))aq and define S x = (f> x /N x . The 
dependence on p(l) has been included in the term we treat as a free normal- 
ization parameter N x , as p{l) is very uncertain both from the theoretical and 
the observational points of view. Although there is a some spread in the pre- 
dictions for the 7-ray flux when coming to specific WIMP models, its spectral 
features are rather generic. As most photons are produced in the hadroniza- 
tion and decay of 7r°s, the shape of the photon spectrum is always peaked 
at half the mass of the pion, about 70 MeV, and it is symmetric around it 
on a logarithmic scale (sometimes this feature is called the u ir° bump", see, 
e.g., [18]). The same is true for the background, but still it may be possible 
to tell signal from background: the signal arises in processes which have all 
the same energy scale, i.e. 2M X , therefore the WIMP induced flux, contrary 
to the background, is spectral index free and shows a sharp cutoff when E n 
approaches the WIMP mass. This is shown in the left panel of Fig. 1, where we 
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plot the differential photon yield per annihilation times the inverse of WIMP 
mass squared, for a few values of the WIMP mass, and assuming that the 
have a single dominant annihilation channel (bb in the case displayed). In the 
same figure, for comparison, the spectral shape of the background is shown: 
as it can be seen, one may hope to identify the WIMP induced component 
as a distortion of the background spectrum at relatively high energies. For a 
given WIMP mass, the photon yields in the different annihilation channels 
are analogous, as shown in the right panel of Fig. 1: solid curves indicate the 
total photon yield, while dashed curves indicate the photon yield in radiative 
processes, i.e. in all processes rather than ir° decays. The spectrum for the it 
and W + W~ channels are very close to the one for bb (differences are mainly 
given by prompt decays before hadronization); only in the t + t~ case, radia- 
tive photon emission is dominant, still with a large bump due to the hadronic 
decay modes of r leptons. 



2.3 EGRET data fit 



EGRET has performed measures in the energy range 30 MeV 10 GeV with 
few bins in the high energy region. Given the paucity of the data in the 
highest end of the energy region in a first approximation it is not sensible 
to discriminate different annihilation channels leading to photons in the final 
state. 

It is then convenient to keep the discussion as general as possible and consider 
a simplified context (a toy-model) with a thermal relic \ f° r which only one 
intermediate annihilation channel is open (Bf = 1 in that channel). Further- 
more we assume the total annihilation cross section to scale with the inverse 
of the relic abundance £l x [26,27]: 

3 • 10 _27 cm 3 s- 1 9fi „ , 

o~v ~ (av) — — 3 • lO-^cmV 1 , (8) 



where (av) is the thermally averaged annihilation cross section. 

Eq.(8) is not valid in presence of resonances or thresholds near the kinemat- 
ically released energy in the annihilation 2M X and of coannihilation effects. 
In the presence of these conditions we can have large deviations from this 
approximate scaling. 

Moreover there are cases in which the inverse proportionality between Q x 
and the annihilation rate only gives a lower bound for the latter since a non- 
thermal relic component can provide the relic density needed to account for 
&CDM [22] . 
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Fig. 2. Fit of the EGRET GC 7-ray data for two sample WIMP models. We fix the 
WIMP mass (M x = 50 GeV in the upper panel, M x = 80.3 GeV in the lower panel) 
and select a single annihilation channel in each of the two cases (bb in the upper 
panel, W~W + in the lower one). Signal and background components are indicated 
separately, while their sum is shown with a solid line. For both models the value of 
the reduced statistical \ 2 variable obtained from the fit is around 5. 

In our toy-model the WIMP mass is kept as a free parameter , as we have 
shown that the photon spectrum is rather sensitive to it. The results we dis- 
cuss below are then dependent on a mass scale and an overall normalization 
parameter and can be easily rescaled to fit an explicit model for which M x 
and av are specified. 

We then fit the EGRET data on the GC excess taken from Table 1 using 
our toy-model. We are not using the two lowest energy bins as for E <C 
1 GeV the background is most probably dominated by inverse Compton and 
bremsstrahlung rather than by 7r° production as we assumed. We repeat this fit 
for different values of the WIMP mass M x and for a few tree-level annihilation 
channels. JV& and N x in Eq.(l) are treated as free parameters which are varied 
in order to minimize the statistical \ 2 variable. The allowed range of variation 
for the background normalization, jVj,, is between 3.2 • 10 20 and 1.8 • 10 21 . The 
two extrema of this variation interval are taken to correspond to a best fit of 
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Fig. 3. Reduced x 2 corresponding to the best fits of the EGRET excess for a WIMP 
model with fixed mass M x and a single annihilation channel allowed (bb is shown 
with a dotted line, W~W + with a solid line); the x 2 has been minimized over the 
normalizations of both the signal and the background (with restrictions as explained 
in the text). Reduced x 2 values should be compared with the value obtained in case 
data are fitted with a background component only, marked with a horizontal line 
in the upper part of the figure. 

the background in agreement with the standard scenario (column 2 in Table 1) 
and to the best fit for the EGRET data from the GC (column 3 in Table 1) 
with N x = 0. 

In Fig. 2 we show our best fits for two sample values of M x and two intermediate 
channels, bb and W + W~ . On the qualitative side, the agreement with the data 
is rather good, even if the reduced x 2 (f° r 6 degrees of freedom) in the two 
examples displayed is still rather large (of the order of 5). This may depend 
on an underestimate of error bars or also on the fact that we are neglecting 
uncertainties in the theoretical predictions for the spectral shapes. It is clear, 
on the other hand, that adding a component due to WIMP annihilations on 
top of the background component greatly improves the agreement between the 
expected spectral form and the one found in the EGRET measurement. This is 
shown in Fig. 3 where the reduced x 2 f° r our Des t fits is shown as a function of 
WIMP mass. These values should be compared with those obtained when only 
the background contribution is included (N x = 0). This case is represented 
by the horizontal line in the figure and gives a reduced x 2 ~ 150. Fig. 3 
also indicates that light WIMP masses are marginally favored over heavier 
WIMP's. Results for other tree-level annihilation states are analogous and 
show the same trends. 



The best fit curves displayed correspond to rather large values of the normal- 
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Fig. 4. Lines of constant reduced \ 2 corresponding to best fits of the EGRET GC 
excess, in the plane WIMP mass M x versus normalization of the WIMP-induced 
signal N x . The plot applies to our toy-model with a single annihilation channel 
allowed, i.e. bb in the case displayed. 

ization parameter N x , a few times 10 4 for the two cases shown in Fig. 3. N x 
tends to increase going to heavier WIMP's as a function inversely proportional 
to the WIMP mass raised to a certain power. This power is a little smaller 
than the value 2 one would have inferred from the fact that the WIMP pair 
density decreases as M~ 2 . An investigation of the dependence of the statistical 
variable x 2 from the WIMP mass and N x is shown in Fig. 4, for the bb annihi- 
lation channel. Given that x 2 = X 2 {Nb, N x , M x ), we minimize it with respect 
to the parameter iVj, for fixed values of N x and M x . The iso level curves for the 
reduced x 2 are plotted in Fig. 4. Fixing a value for the WIMP mass and going 
from the left-hand side to the right-hand side of the figure, we move from the 
case in which the WIMP signal is marginal with respect of the background 
flux to the case where the sum of the two reproduces most closely the data. 
Further increasing the value of N x we reach a region where the WIMP signal 
exceeds the flux detected by EGRET. 

We recall that in our toy-model N x is identified with the halo model depen- 
dent function (J('0))aq f° r the Galactic center direction ip = 0° and EGRET 
angular acceptance AQ = 10~ 3 sr. As already mentioned, the distribution of 
dark matter in the inner part of the Galaxy is still a controversial issue. Dy- 
namical measurements show that dark matter is, in the Galactic center region, 
just a subdominant component with respect to baryonic matter, but lack the 
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Table 2 

Values of {J(0)) An for two different Ail's and for three different density profiles, 
see the text for details. 



Profile (J(0)) A n (AO = KT 3 sr) <J(0)) An (Afi = lO" 5 sr) 



Navarro, Frenk, White 


1.21 


10 3 


1.26 • 


10 4 


Moore et al. 


1.05 


■ 10 5 


9.46 • 


10 6 


Modified isothermal 


3.03 


10 1 


3.03 • 


10 1 



resolution we need for an estimate of (J). On the other hand, N-body sim- 
ulation of structure formation in a CDM Universe, find dark matter density 
profiles which are singular toward the center of the Galaxy, possibly scaling 
as 1/r (profile of Navarro, Frenk & White [28] (NFW)) or l/r L5 (Moore et 
al. profile [29]) as the galactocentric distance r — > 0. The corresponding values 
of (J), and hence of N x are listed in the second column of Table 2 assuming, as 
commonly done, that the dark matter density at the Sun's galactocentric dis- 
tance is equal to 0.3 GeVcm -3 . One should however keep in mind that it has 
been questioned whether the NFW and the Moore et al. profiles can be used 
to describe inner dark matter halos (especially for smaller galaxies, see, e.g., 
[30] for a review). These halo models give a snapshot of the Galaxy before the 
baryon infall. The appearance of a massive black hole at the Galactic center 
and of the stellar components may sensibly modify such pictures with further 
enhancements (but a depletion is possible as well) of the central dark matter 
density [31,32]. For comparison, in Table 2 we give the value of (J) for the 
modified isothermal sphere profile, which is non-singular toward the Galactic 
center and, as well known, give a normalization for the dark matter induced 
fluxes well below the background and the sensitivity of even next-generation 
detectors. All three halo profiles listed in the Table are consistent with avail- 
able dynamical constraints on the Galaxy. We conclude then that N x can be 
at the level needed in our toy-model to reproduce the EGRET excess, but at 
the same time that both larger or smaller values are feasible as well. 



3 EGRET Excess as Mapped by GLAST 



Much more information on the nature of the EGRET excess at the Galactic 
center will be provided by the GLAST telescope. With respect to EGRET, 
GLAST will cover a wider energy range, with an increased effective area and 
better energy and angular resolution. Besides pinning down features in the 
WIMP-induced flux mediated by n° decays, GLAST will have the power to 
search for the monochromatic gamma-ray flux eventually arising from the pair 
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Fig. 5. Simulation of the data set which will be obtained with GLAST in 2 years, 
in case the EGRET GC excess is due to the WIMP-induced flux shown in one of 
the sample fits in Fig. 2 (lower panel). The error bars refer to statistical errors for 
the chosen energy binning and for the angular acceptance AQ = 10 -3 sr. 

annihilation, at 1-loop level, of non-relativistic WIMP's into a two-body final 
state containing a photon (for neutralinos as WIMP dark matter candidates, 
two such states are allowed: XX ~^ %1 [33] and XX ~^ 1% [34], producing 
photons with energy equal to, respectively, M x and M x (1 — Mf/4M 2 )). The 
discussion of GLAST potential to detect the monochromatic components is 
postponed to a future analysis; we focus here on the term with continuum 
energy spectrum. 

We start by supposing that the GC excess as mapped by EGRET is indeed 
due to WIMP annihilations and extrapolate what kind of data GLAST would 
collect about it. For the performance of the detector, we rely on a simplified 
picture emerging from the latest simulations [14]. We assume that, on average 
in the energy interval of interest to us, the instrument has an energy resolution 
of 10%, angular resolution of 1CT 5 sr (~0.1°), and a peak effective area of 
11000 cm 2 . GLAST will perform an all sky survey, rather than operating in 
the pointing mode. We assume a data acquisition time of 2 years and derive 
the fraction of time the GC center is visible by the instrument. We simulate a 
sky survey with ± 35° rocking and take into account the loss of exposure due 
to the South Atlantic Anomaly (SAA) passages, which corresponds to about 
14.2% of the orbiting time. We find that the fraction of time that the GC is 
within the GLAST field of view is equal to 0.592, and that the net fraction 
of time that the source can be observed is 0.508; the reduction in effective 
area for sources which are not located at the instrument zenith gives a mean 
effective area equal to 60% of the peak effective area. 
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Fig. 5 shows a projection for the GC flux which could be measured by GLAST, 
assuming the spectrum and normalization for the dark matter source and the 
normalization for the background as derived from the fit of the EGRET data 
in Fig. 2 (lower panel). The simulated data points are derived by assuming 
as angular acceptance the EGRET angular resolution (Af2 ~ 1CT 3 sr), and 
choosing the energy bin widths to be of the order of 10% of their central values 
to take advantage of the GLAST energy resolution. The error bars displayed 
are associated to the statistical error only. If we try to fit the simulated data 
with the spectral shape of the background and a free normalization, the re- 
duced x 2 we get is higher than 2 • 10 3 , much larger than the reduced x 2 we 
obtained in the corresponding fit for the EGRET data set. 

GLAST will also allow to search for an eventual angular signature of the dark 
matter source. A 10" 5 sr angular resolution implies that the telescope will 
map the GC resolving regions with a precision of 7 pc (assuming the sun 
galactocentric distance is 8.5 kpc). 

Most models for the distribution of dark matter in the Galaxy such as the 
NFW and the Moore et al. models, predict an enhancement in the dark matter 
distribution toward the Galactic center on a scale which exceeds this size. 
There is then the chance that WIMP annihilations in the GC region may 
produce a detectable flux on an angular scale exceeding the angular resolution 
of GLAST. It is then sensible to investigate whether GLAST will detect the 
flux in the example shown in the Fig. 5 as coming from a point source located 
at the GC or from a diffuse source with degrading intensity increasing the 
angle between the line of sight and the GC direction. To do that we need 
to focus on a specific model for the dark matter density profile in the inner 
Galaxy. Inspired by the NFW and the Moore et al. profiles, we assume that 
the WIMP density toward the GC scales like: 



We normalize p(r) by fixing the local WIMP density, i.e. the density at the 
galactocentric distance r = 8.5 kpc, to be equal to p = 0.3 GeV/cm 3 . 7 is 
kept as a free parameter. To avoid the singularity in r = we have introduced 
a lower cut-off r m j„ = 10 -5 kpc, corresponding to a distance from the GC 
below which we assume that the power law behavior cannot be trusted. 

For the sample toy- model shown in Fig. 2 (lower panel) and in Fig. 5, N x = 
(J(0)) A n(Afi = 10 -3 sr) = 8.5 • 10 4 from which we infer that 7 = 1.54. Values 
of (J(i()))aci, for such 7 and for the dark matter density profile specified in 
Eq.(9), are shown in Fig. 6 as a function of ip and for a few values of Afl. 
We can then calculate the expected flux for the model obtained in the fit of 




r > r. 



r < r. 



mm 



(9) 
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Fig. 6. Angular dependence for the WIMP signal displayed in the lower panel of 
Fig. 2 and in Fig. 5 in the case in which the dark matter density profile p(r) has 
the power law form introduced in Eq.(9). ip is the angle between the direction of 
observation and that of the GC. (J(Y>)}An coincides with N x for the toy-model we 
introduced. Three sample angular acceptances AQ are considered. 




Fig. 7. For the WIMP-induced flux corresponding to the model shown in Fig. 5 and 
Fig. 6, we compute the expected flux that would be obtained by GLAST in 2 years. 
We fix the angular acceptance AO and denote the angle between the direction of 
observation and the GC by tp. In the figure we plot the reduced \ 2 obtained by 
fitting the expected flux with a background component only. We find that x 2 2> 1 
for ip much larger than the angular resolution of GLAST (vertical dotted line in the 
figure). We then conclude that GLAST will be able to resolve the angular structure 
of the signal for the case considered. 



15 



EGRET data, fixing the direction of observation and the angular acceptance. 
Finally we compute the expected flux that, with the above provisions, GLAST 
will collect. We have also included a background component independent from 
tp which is equal to that we have estimated from the fit of the data by EGRET. 

In Fig. 7, as a function of the direction ip and for a few values of AQ, we 
plot the value of the reduced x 2 we get by fitting in each case the expected 
flux with a component that has the spectral shape of the background and free 
normalization. As x 2 ^ 1 on angular scales much larger than the angular 
resolution, this indicates that GLAST will resolve the angular structure of the 
dark matter source at the GC. 

iFrom the figure we see also that, for ip — 0, the largest x 2 are obtained for 
the minimum Afl considered, i.e. an angular acceptance equal to the angular 
resolution of the detector. This indicates that, for the specific halo model 
considered, the ratio of the signal to square root of the background increases 
going to smaller and smaller angular acceptances. 

Values of the reduced x 2 have been obtained so far by supposing that the spec- 
tral shape of the background is known. Actually, what we have done is to fix it 
according to one of the currently favored scenarios. Slight discrepancies with 
respect to this model are plausible. On the other hand, GLAST will perform 
an all sky survey in which the background component will be accurately mea- 
sured at all longitudes and latitudes. It might still be problematic to choose 
the background normalization in the Galactic center direction, if an excess is 
indeed found in that direction. But it will be possible to relax the assumption 
that the spectral shape of the background is known from theory, as it will be 
possible to extrapolate it from the data at higher latitudes and longitudes. 
A (measured) spectral shape which is different from what we assumed would 
slightly change the numerical predictions we derived so far; on the other hand 
the general features would remain exactly the same. Keeping this in mind, one 
should not take a reduced x 2 of 1 as a strict discriminator, e.g. in Fig. 7, to 
decide whether the signal could be resolved from the background. 



4 GLAST Performance for a Weaker Source 

Although we have shown that the flux measured by EGRET is compatible with 
being due to WIMP annihilations, until more accurate data are available it will 
not be possible to discriminate this solution from other plausible scenarios. In 
particular, given the rather poor angular resolution of the EGRET instrument, 
there is even the possibility that the EGRET excess is not actually associated 
to a source located at the Galactic center. This is the case, e.g., if the flux 
is identified with inverse Compton emission from the electrons responsible for 
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the synchrotron emission in the radio arc around the Galactic center [35]. 
Should the next generation of gamma-ray observations solve the puzzle in this 
direction, there still would be room for finding a component due to WIMP 
annihilations in the Galactic center region, maybe associated to a dark matter 
source weaker than the one we postulated so far. Hence we extend the analysis 
we performed to investigate the potential of GLAST to single out such a 
source. 

As we have shown, for a halo profile of the type given in Eq.(9), it is advanta- 
geous to focus on a region which is as small as possible around the GC. Hence 
we consider a survey with angular acceptance equal to the angular resolution 
of GLAST, ~ 10~ 5 sr. We assume that the background component is still due 
to diffuse emission from primary cosmic rays interacting with the interstellar 
medium. Hence we keep the spectral form implemented so far for Sb, with 
a normalization such that it matches at least the higher latitude measured 
flux, i.e. Nb > 3.2 ■ 10 20 for a background at least at the level of the flux 
reported in the second column in Table 1. For a given WIMP model, i.e. fix- 
ing S x , we search then for the minimum ratio between the two normalization 
factors N x /N b that is needed to eventually discriminate with the GLAST tele- 
scope the WIMP annihilation signal from the background. N x is now equal 
to (J(0))aq(AJ1 = 10~ 5 sr): sample values for this quantity are listed in the 
third column of Table 2 for the three halo profiles introduced in Section 2.3. 
Again, the huge spread in the predictions (possibly even further amplified by 
the redistribution of dark matter particles during the formation of the cen- 
tral black hole [31,32]) reflects our lack of knowledge about the dark matter 
density in the Galactic center region. 

For each pair N x and Nb, assuming the same instrumental performance, ex- 
posure time and energy binning specified in Section 3, we simulate the corre- 
sponding data set that GLAST would obtain, i.e. the expected flux measure- 
ments with the associated statistical error for the chosen energy binning, anal- 
ogously to what is shown in Fig. 5 for the EGRET GC source. The criterion 
to discriminate whether the WIMP component would be singled out is based 
on the usual x 2 test statistic. We have computed the reduced x 2 between the 
number of counts expected in each energy bin for the two hypothesis: WIMP 
signal plus background and background only. Taking into account the number 
of degrees of freedom, which in our case is equal to the number of energy 
bins, the signal plus background curve is distinguishable from the background 
only curve, for a reduced x 2 > constant. This constant in uniquely deter- 
mined by the number of degrees of freedom and by the confidence level we 
want to reach. We have also checked our results against those obtained with 
the likelihood ratio method [36,37], obtaining no discrepancies. 2 This latter 



2 We thank G.Ganis for providing the software package we used to perform this 
analysis. 
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Fig. 8. Minimum ratio between the normalization of the WIMP signal N x and the 
background normalization Nf, such that the WIMP induced signal would be singled 
out of the background with GLAST. We are referring to a toy-model with a single 
annihilation channel allowed, i.e. bb in the case displayed. In the linear regression 
fit marked in the figure y = log 10 (jf-j and x = log 10 M x . 




method is especially suited for the case we have at hand: to decide whether 
a certain event belongs to the background only hypothesis (Hq) or to signal 
plus background hypothesis (Hi), one starts by constructing two probability 
distributions, P and Pi, for an estimator F = L(Hi)/L(H ), which is the 
ratio between the likelihoods L of the two hypotheses. In our case, since we 
are interested in counting, we can choose the Poisson distribution to obtain 
the likelihood. Comparing the two distributions one can decide, at a certain 
confidence level, if they will result distinguishable or not, once it is fixed the 
accuracy of the experimental data that will be used for the discrimination. 
The likelihood ratio method is in general more powerful than the x 2 one > 
since, in addition to giving the probability of a certain set of data to belong to 
the signal plus background probability distribution, it allows to compute the 
probability to be wrong when accepting such hypothesis, the so called power 
of the test, considering the background only hypothesis as the true one. 

As a sample test case, we consider again the toy-model introduced in Sec- 
tion 2.3 for a single WIMP annihilation channel, e.g., bb. In Fig. 8 we plot, as 
a function of the WIMP mass, the minimum value of the ratio N x /Nb needed 
for a discrimination of the WIMP signal from the background. The corre- 
sponding values of N x = (J(0))aq(AJ1 = 1CT 5 sr) are of the order of ~ 10 3 
for a WIMP mass of about 50 GeV. They increase approximately as M* -8 for 
heavier particles. These values are larger than the ones one would obtain from 
a smooth profile, see Table 2. A local enhancement in the WIMP dark matter 
density at the GC is then required to match our values of N x . Such enhance- 
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merit seems to be smaller than the one needed to fit the EGRET excess in the 
previous section (iV x of the order of ~ 10 4 -~ 10 5 ). As a word of caution we 
remark that the values of N x we obtained in the two analysis should not be 
directly compared in principle as the angular acceptance in the two cases is 
different: Af2 = 1CT 5 sr in the current case, while we had Af2 = 1CT 3 sr when 
we fitted the EGRET data set. A halo model has to be specified to translate 
one into the other. 



5 A Specific WIMP: the Lightest Neutralino in the mSUGRA 
Framework 

As a sample application of the generic tool we discussed so far, we focus now on 
the most widely studied WIMP dark matter candidate, the lightest neutralino, 
in the most restrictive supersymmetric extension of the Standard Model, the 
minimal supergravity (mSUGRA) framework [38]. This setup has been con- 
sidered extensively in the contest of dark matter detection (a list of recent 
references includes, e.g., [39,40,41,42,43,44]) and therefore the comparison of 
our results with previous work and other complementary techniques should 
be transparent in this case. 

In the general framework of the minimal supersymmetric extension of the 
Standard Model (MSSM), the lightest neutralino is the lightest mass eigen- 
state obtained from the superposition of four interaction eigenstates, the su- 
persymmetric partners of the neutral gauge bosons (the bino and the wino) 
and Higgs bosons (two Higgsinos). Its mass, composition and couplings with 
Standard Model particles and other superpartners are a function of the sev- 
eral free parameters one needs to introduce to define such supersymmetric 
extension. In the mSUGRA model, universality at the grand unification scale 
is imposed. With this assumption the number of free parameters is limited to 
five: 

m i/2, m o, sign(fi), A and tan/5, 

where mo is the common scalar mass, m.1/2 is the common gaugino mass 
and Aq is the proportionality factor between the supersymmetry breaking 
trilinear couplings and the Yukawa couplings, tan/5 denotes the ratio of the 
vacuum expectation values of the two neutral components of the SU(2) Higgs 
doublet, while the Higgs mixing \i is determined (up to a sign) by imposing 
the Electro- Weak Symmetry Breaking (EWSB) conditions at the weak scale. 
In this context the MSSM can be regarded as an effective low energy theory. 
The parameters at the weak energy scale are determined by the evolution of 
those at the unification scale, according to the renormalization group equations 
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(RGEs). For this purpose, we have made use of the ISASUGRA RGE package 
in the ISAJET 7.64 software [45]. 

After fixing the five mSUGRA parameters at the unification scale, we extract 
from the ISASUGRA output the weak-scale supersymmetric mass spectrum 
and the relative mixings. Cases in which the lightest neutralino is not the light- 
est supersymmetric particle or there is no radiative EWSB are disregarded. 
The ISASUGRA output is then used as an input in the DarkSUSY package[24]. 
The latter is exploited to: 

• reject models which violate limits recommended by the Particle Data Group 
2002 (PDG) [46]; 

• compute the neutralino relic abundance, with full numerical solution of the 
density evolution equation including resonances, threshold effects and all 
possible coannihilation processes [47]; 

• compute the neutralino annihilation rate at zero temperature in all kine- 
matically allowed tree-level final states (including fermions, gauge bosons 
and Higgs bosons; 

• estimate the induced gamma-ray yield by linking to the results of the sim- 
ulations performed with the Lund Monte Carlo program Pythia [23] as 
implemented in the DarkSUSYpackage. 

Note that none of the approximations implemented in the toy-model intro- 
duced above, regarding the estimate of the relic density or the annihilation 
cross section, are applied here. 

We are ready then to exploit the procedure outlined in Section 4 to study in 
what region of the mSUGRA parameter space, for a given dark matter halo 
profile, the induced continuum 7— ray flux would be detectable by GLAST. 
Fixing tan /3, A and sgn(/i), we have performed a scan in the (m , mi/ 2 ) plane 
searching for the minimum dark matter density, in the GC region, needed to be 
able to single out the neutralino annihilation signal with GLAST. To do this 
we have followed the same discrimination criteria described in Section 4 that 
we recapitulate here for the reader's benefit. First we estimate the statistical 
error (la) on GLAST data to be the square root of the number of events. To 
compute the latter we multiply the flux by the effective area of the detector, 
by the total observational time and the angular resolution AQ = 10~ 5 sr. 
Then for each value of the pair (mo, 7711/2) we compute the difference between 
the fluxes 7 = ft + <f> x = N b Sb + N X S X and ft = ft = N b Sb- If 7 — 7 > 3a 
we consider the SUSY model with those values of (m ,m 1 / 2 ) to be detectable 
by GLAST. In Figs. 9 and 10 we show the isolevel curves for the minimum 
allowed value of N x = (J(0))aq(AQ = 10~ 5 sr) for the signal detection, in the 
(mo,mi/ 2 ) plane and for five sample sets of the other parameters. 

The colored regions in the figures represent portions of the parameter space 
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Fig. 9. Contour plots in the mSUGRA (?n<o> m i/2) plane, for the value of the nor- 
malization factor N x , that allows the detection of the neutralino 7 ray signal, with 
GLAST. In the green region 0.13 < n x h 2 < 0.3, while the red region corresponds to 
the WMAP range 0.09 < Q x h? — O--^ [1]- The black region corresponds to models 
that are excluded either by incorrect EWSB, LEP bounds violations [46] or because 
the neutralino is not the LSP. In the dark shaded region rrih < 114.3 GeV [46], 
where ho is the lightest Higgs. 

where the neutralino has the right cosmological relic density to constitute 
CDM: in the green region 0.13 < Q x h 2 < 0.3, while the red region corresponds 
to the WMAP range 0.09 < VL x h 2 < 0.13 [1]. It is interesting to remark that in 
the cosmo logically favored regions the values of (J) which allow the detection 
of the WIMP signal by GLAST are among the smallest in the interval of 
variability of the variable (J) itself. To some extent this was anticipated in 
Eq. (8) that was one of the approximations of our toy-model: low values of the 
relic abundance lead to high values of the annihilation cross section. 

The cosmologically favored regions correspond to two regimes which have been 
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Fig. 10. Contour plots in the mSUGRA (mo, ^1/2) plane for tan (3 = 55. Left panel: 
values of the normalization factor N x , that allow the detection of the neutralino 7 
ray signal, with GLAST. Right panel: values of the neutralino mass. The excluded 
and colored regions are as in Fig. 9. 

extensively studied in the literature. The first regime is given by models where 
the neutralino tends to be a very pure bino. It comprises the region where the 
neutralino is light (lower parts of each panel, compare with Fig. 12) and the 
region where the stau coannihilation is active (the low mo region on the left of 
the first three panels, see, e.g., [48]). The second regime, sometimes dubbed as 
"focus-point" region [49,50], is the region on the right-hand side of the panels 
with tan (3 equal to 10, 50 and 55, close to the region in parameter space where 
there is no EWSB, in which the neutralino has a relevant Higgsino component. 

Comparing the values of (J(0))An(Aft = 10~ 5 sr) obtained with those listed in 
Table 2, we see that in the case of a moderate enhancement of the dark matter 
density toward the GC as predicted by the NFW profile or an even slightly 
milder singular profile, there is a fair portion of the mSUGRA parameter space 
which gives fluxes detectable with GLAST. It is worth to observe that high 
tan/3 models (i.e. tan/3 = 50 and tan/3 = 55) possess the largest cosmologi- 
cally favored regions. Furthermore they give the highest chances to single out 
the 7-ray signal from neutralino annihilations with GLAST. We can compare 
our results with those of [39]. There it is assumed that GLAST has a certain 
sensitivity to the integrated continuum 7— ray flux from a region around the 
GC, of an extension 10 2 times wider than the GLAST angular resolution. The 
neutralino-induced signal is then supposed to be detectable if its integrated 
flux exceeds such sensitivity. In Figs. 18 and 19 of [39] are given the regions 
which can be probed by GLAST in case of ( J(0)) An (Aft = 10" 3 sr) = 500. 
They are in qualitative agreement with our corresponding predictions (first 
and fourth panels of Fig. 9) for Aft = 10~ 5 sr. 
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Fig. 11. Contour plots in the mSUGRA (w-0i m i/2) plane, for values of N x that 
are already excluded by EGRET data at 5a confidence level. In the green re- 
gion 0.13 < Q x h 2 < 0.3, while the red region corresponds to the WMAP range 
0.09 < Q x /i 2 < 0.13 [1]. The black region corresponds to models that are excluded 
either by incorrect EWSB, LEP bounds violations [46] or because the neutralino 
is not the LSP. In the dark shaded region rrih Q < 114.3 GeV [46], where ho is the 
lightest Higgs. 

In addition to this study of the GLAST sensitivity, we have tried to single out 
the regions of the mSUGRA parameter space ((m ,m 1 / 2 ) for fixed tan(/3), A 
and sign(n)), which are already experimentally excluded, due to a neutralino- 
induced 7-ray flux exceeding the GC EGRET data of table 1. In Fig. 11 we 
show isolevel curves for the maximal (J(0))aq(AH = 10~ 3 sr) marginally 
consistent with the data at the 5a level. We can observe that the values we 
obtain are generally significantly larger than the corresponding sample values 
in Table 9. This result implies that rather weak constraints on the mSUGRA 
parameter space can actually be imposed, on the basis of the GC 7— ray flux 
measured by EGRET. 
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Fig. 12. Neutralino mass contour plots in the (mo,mi/ 2 ) mSUGRA plane. The 
black region corresponds to models that are excluded either by incorrect EWSB, 
LEP bounds violations [46] or because the neutralino is not the LSP. In the dark 
shaded region < 114.3 GeV [46], where ho is the lightest Higgs. 



One could finally wonder how the GC EGRET data could be fitted in the con- 
text of the mSUGRA models. In order to answer to this question, we rely on 
the general analysis that we have performed in the toy-model scheme (see sec- 
tion 2). In particular we could extract from Fig. 4 the values of the neutralino 
mass that allow the best fit of the GC EGRET data 3 . Looking at Fig. 12 we 
can then single out the (m , ^1/2) parameter regions that corresponds to such 
values. 



3 Recall that the possibility to fit EGRET data is quite insensitive to the dominant 
annihilation channel. 



24 



6 Conclusions 



We have found that the excess in the 7-ray flux detected by the EGRET tele- 
scope toward the Galactic center shows spectral features which are compatible 
with an exotic component due to WIMP annihilations, especially for WIMP 
masses in the lower end of the mass range currently considered for WIMP dark 
matter candidates. For the WIMP-induced flux to be at the level of the mea- 
sured flux, a fairly large dark matter density is needed in the Galactic center 
region; indeed, such density enhancements are found in N-body simulations of 
halo profiles in cold dark matter cosmologies. 

Although it is not possible with present data on the Galactic center excess 
to discriminate between the interpretation we propose here and other viable 
explanations, we have shown that, with the data that will be collected by the 
GLAST, the next major 7-ray telescope in space, it will be possible to identify 
both spectral and angular signature expected for a WIMP-induced component. 
If on the other hand the data will point to an alternative explanation, there 
will still be the chance for the GLAST telescope to single out a (weaker) 
dark matter source. The potentials of GLAST have been explored both in 
the contest of a generic simplified toy-model for WIMP dark matter and for 
one of the most widely studied WIMP dark matter candidate, the lightest 
neutralino, in the minimal supergravity framework. We find that even in case 
of moderately singular dark matter profiles, there are regions in the parameter 
space which will be probed by GLAST, especially in the high tan (3 case. We 
find, on the contrary, that limits from current EGRET data are rather weak. 
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